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Abstract 

We reprove the A-Lemma for finite dimensional gradient flows by gen¬ 
eralizing the well-known contraction method proof of the local (un)stable 
manifold theorem. This only relies on the forward Cauchy problem. We 
obtain a rather quantitative description of (un)stable foliations which al¬ 
lows to equip each leaf with a copy of the flow on the central leaf - the 
local (un)stable manifold. These dynamical thickenings are key tools in 
our recent work [Web]. The present paper provides their construction. 
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1 Introduction and main results 

Throughout (M",g) denotes a Riemannian manifold of finite dimension n and 
/ is a function on M of class (7’’+^ with r > 1. We assume further that the 
downward gradient equation 

Ht) ■= f) ° u(t) ( 1 ) 

for curves u : [0, oo) —)■ M generates a complete forward flow This 

holds true, for instance, whenever / is exhaustive, that is whenever all sublevel 
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Figure 1: Local coordinates near non-degenerate critical point 


sets M“ = {/ < a} are compact. By x we denote throughout a non-degenerate 
critical point of / of Morse index k = ind(a;) with 1 < fc < n — 1; for A: S {0, n} 
there is no A-Lemma. In the following we give an overview of the contents of 
this paper. Definitions and proofs are given in subsequent sections. 

The A-Lemma^ was proved by Palis [Pal67,Pal69] in the late 60’s. Its back¬ 
ward version asserts that, given a hyperbolic singularity x of a vector field 
X on a finite dimensional manifold, the corresponding backward flow applied 
to any disk D transverse to the unstable manifold kF“(x) converges in and 
locally near x to the stable manifold of x; cf. [PdM82, Ch. 2 §7] and Figure 2. 
Since the A-Lemma is a local result we pick convenient coordinates about x. 

Definition 1.1 (Local coordinates and notation - Figure 1). 

(HI) Associated to the Hessian operator A = D{yf)^ there is the spectral 
splitting X := T^M = E~ © in (7) with spectral projections 7r± and 
where E~l^ = T^W^I’^ix). Consider the spectral gap d > 0 of A, see (6), 
and fix a constant A G (0,(i) - which will become the rate of exponential 
decay. Consider the constants 0 < d < /i < d defined by Figure 4. 

(H2) Express the downward gradient equation (1) via exponential coordinates 
in the form of the ODE (12) on a ball Bp^ C X := T^M whose radius 
Po < 1 satisfies (14). In these coordinates the origin represents x and the 
local flow will be denoted hy (j) = The non-linearity h of this ODE, 
also shown in (12), satisfies the Lipschitz estimates of Lemma 2.5. 

Notation. The notation of coordinate representatives of global objects 
such as W‘‘{x) C M will be the global notation without the x such as IT®. 

(H3) By Remark 3.9 we take our local model of the form V = D x where 
V C Bpg and C represent a/de/scending disks, see (18) 

and Remark 3.6, whose parameter 2c is determined by the inclusions (25). 
In this local model the graph maps (20) and (23) of the local un/stable 
manifolds take on the rather simple form 

7"“ = (id, 0):D- ^ D- X D+, = (0, id) : D+ ->■ D- x D+. 

^ The naming inclination or A-Lemma is due to the fact that it asserts convergence 
of a family of disks to a given one in the (un)stable manifold, so the inclination of suitable 
tangent vectors converges. In [Pal69] these inclinations were denoted by A^. 
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Figure 2: Backward A-Lemma 


(H4) Suppose £ S (0,^) and consider the descending sphere S'" = dW^ con¬ 
tained in = D~. Pick a constant x S (0,1] such that the radius x 
ball 5+ C E'^ lies in the ascending disk C = iA+. Consequently 
the hypersurface I? := S" x is contained in the interior of V. The 
fiber Va := {a} x over a € S” is a translated ball of codimension k. 

Theorem 1.2 (Backward A-Lemma). Given the local model on V = x 
C E~ 0 E'^ = T^M provided by Definition 1.1 nearby a non-degenerate 
critical point x of a function / : M —)■ K, pick e G (0,<^). Then the following 
is true. There is a ball B~^ C E~^ about the origin of some radius ^ > 0, a 
constant Tq > 0, and a Lipschitz continuous map 

g : [To, oo) X S" X B+ ^ y C T- © E+ 

{T,z_,z+)^ (g'^_{z+),z+^ =-.gJ_{z+) 

defined by (36) below which is of class C" in z_ and z+. It satisfies the identity 
= 4’-t(z-) and the graph of Gf_ consists of those z gV which satisfy 
TT+z G and reach the fiber = {z-} x at time T, that is 

GL iB+) = n (E- X B+) . 

Moreover, the graph map gj_ converges uniformly in , as T -G oo, to the 
local stable manifold graph map g°° as illustrated by Figure 2. More precisely, 

\\dgj_{z+)v - dg°°{z+)v\\ = ||dGf_ ( 0 +)z;|| < ||n|| , 

for all T > Tq, z- G Sf, z+ G B'^, and v G T+ where c* > 0 is a constant. 
Estimate two requires f G G^^ near x in which case T -^g{T, z_, z+) is also 
Lipschitz continuous.'^ 


^ The condition / E near x (satisfied for r >2) hinges on the Lipschitz Lemma 2.5. 
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Remark 1.3 (Contraction method proof of A-Lemma). The contraction method 
proof presented below has its own dis/advantages. On the worrying side, for 
/ G that is X = —Vf G we only obtain C® convergence so far, whereas 
we do get convergence if / is of class near x, e.g. if / G On the 
bright side, we get rather useful quantitative control on each of the involved 
variables such as time t, the variable describing the dislocated disks, and their 
dependence on the base point; for details see Theorem 1.2. Most importantly, 
this quantitative control lends itself to construct foliations around critical points, 
foliations associated to global objects, namely sublevel sets, and to equip each 
leaf with its own semi-flow; see [Webl4b,Web] for first applications. 

Remark 1.4. Theorem 1.2 is a family backward A-Lemma, namely for the fam¬ 
ily of disks T> = X Given a general hypersurface V whose intersection 
with the unstable manifold is transverse and equal to S'", one needs to change 
coordinates to bring into the required normal form 5+. To achieve this 
apply the parametrized Implicit Function Theorem to construct a C'' diffeomor- 
phism of V which is the identity outside a small neighborhood U of the compact 
intersection locus S'" and maps the part of V near S" to S“ x 5+ for some x. 

We substituted the backward flow of the disk T>z_ by its pre-image 4>t~^'Dz_ , 
because there are interesting situations, cf. [HenSl], with only a forward semi¬ 
flow. For instance, the heat flow on the loop space of a closed Riemannian 
manifold [SW06] is just a forward semf-flow. Surprisingly, although there is no 
backward flow, there is a backward A-Lemma [Webl4a] under which the hyper¬ 
surface T) moves backward in time and foliates an open set. Furthermore, it is 
even still possible to construct a non-trivial Morse complex [Web 13, Web 14b]. 
Consequently the - in finite dimensions - common situation of a genuine for¬ 
ward and backward flow with stable and unstable manifolds, all being finite 
dimensional, allows for (too) many choices. In order to emphasize the necessary 
elements and, at the same time, to introduce the reader tacitly to the infinite 
dimensional heat flow scenario, we subject ourselves to the following convention. 

Convention 1.5 (No backward Cauchy problem). We shall use existence of a 
solution to the Cauchy problem only in forward time and ignore the back¬ 
ward Cauchy problem alltogether. However, on the unstable manifold there is 
still a natural procedure to define a backward flow; see Definition 3.3. Since 
there is no Cauchy problem involved, we will (and need to) use this so-called 
algebraic backward flow; this is coherent with the infinite dimensional case 
cited above. In finite dimensions the algebraic backward flow coincides with the 
one obtained via the Cauchy problem. So we do use a backward flow, but with¬ 
out solving the backward Cauchy problem and only along the unstable manifold. 
To put things positively, our self-imposed lack of mathematical structure elim¬ 
inates the possibility of choices and therefore reveals those elements which are 
essential to define a Morse complex - or even just an unstable manifold. 

The simple but far reaching idea, see [Web 14a, Web 14b], which avoids any 
backward Cauchy problem is to look at the pre-image under the time-T-map of 
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Figure 3: Conley pair (7V®,L®) associated to critical point x 


a given disk ^^{x) transverse at q to the unstable manifold W'^{x). The obvious 
but non-trivial problem is then to show that these sets are in fact submanifolds 
of M. To show this we shall write them locally near x as graphs over the 
stable tangent space E~^ = TxW‘^{x) using the contraction method which 
in this case consists of interpreting the backward Cauchy problem as a “mixed 
Cauchy problem”. The latter can be reformulated in terms of the fixed point of 
a contraction, both depending on parameters. Thereby we obtain a new proof 
of the A-Lemma in finite dimensions; see [Webl4a] for the infinite dimensional 
semi-flow case of a parabolic PDF where this program has been carried out first. 
A key point of interest will be to control how the pre-images depend on time T. 

Stable foliations 

Borrowing from [Sal90] we define for each choice of constants e, r > 0 a pair 

K = N^{e,T) : = {pe M \ f{p) <c + e, f{prP) >c-£}^, (3) 

where {.. .}x denotes the path connected component that contains cc, and 

L'l = LI{£, t) := {p e iV® I f{p 2 rp) <c-£}. (4) 

Figure 3 shows a typical such pair, illustrates what is called the exit set property 
of L®, and indicates hypersurfaces which are characterized by the fact that each 
point reaches the level set {/ = c — e} in the same time. The points on the 
stable manifold never reach level c — e, so they are assigned the time label 
00 . By the backward A-Lemma, Theorem 1.2, these hypersurfaces foliate some 
neighborhood of x. However, the neighborhood and so the leaves have no global 
meaning so far. It is the content of Theorem 5.1 that the pairs {Nf., L%) which 
are defined in terms of level sets are foliations; for a first application see [Web]. 

Remark 1.6 (Forward A-Lemma and unstable foliations). Since the backward 
A-Lemma as well as the stable foliations (iV®, L®) with respect to / are of local 
nature at x one can simply obtain the forward version of Theorem 1.2 and the 
unstable foliations {N^,Lf.) with respect to / by taking the former ones with 
respect to —/. (If necessary, cut off —/ to obtain a complete forward flow.) 
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2 Preliminaries 


Throughout we denote by V = V® the Levi-Civita connection associated to the 
Riemannian metric g on M which we also denote by (•,•). Given a singularity 
p of a vector field X on M we denote by DXp linearization dXp : TpM —)■ 
Tx(p)TM = TpM(BTpM followed by projection onto the second component. The 
isomorphism is natural, because Xij}) = 0. There is the identity DXp = {X.X)p. 

Hessian and spectrum 

By Grit/ we denote the set of critical point of /. Given p € Grit/, in other 
words a singularity oi X = V/, then there is a linear symmetric operator 
A = Ap : TpM —)■ TpM uniquely determined by the identity 

Hessp/(u, ru) = {Av,w) (5) 

for all v,w G TpM. The bilinear form HesSp/ is defined in local coordinates in 
terms of the (symmetric) matrix with entries {d^ f /dx^dx^){p). It holds that 
Av = D{Xf )pV = iXyXf)p where v is a vector field on M with Vp = v. By 
symmetry the spectrum of the Hessian operator A is real. If it does not 
contain zero, then p is called a non-degenerate critical point of /. 

Suppose a; is a non-degenerate critical point of /. Then by the spectral 
theorem for symmetric operators the spectrum of ^ = A^ is given by n reals 

Ai < • • • < Afc < 0 < Afc+i < • • • < A„ (6) 

counting multiplicities. The distance d > 0 between 0 and the spectrum of 
A is called the spectral gap of A. The number k of negative eigenvalues is 
called the Morse index of the critical point x and denoted by ind(a;). Denote 
by E~ and E~^ the span of all eigenvectors associated to negative and positive 
eigenvalues, respectively. Then there is the orthogonal splitting 

T^M = E- ®E+ (7) 

with associated orthogonal projections 7r_ : T^M —)■ E~ and 7r+ : T^M —)■ 
E'^ . Since A preserves E~ and E'^ it restricts to linear operators A~ := ^o7r_ 
and A+ := Aoxj^ on E~ and respectively. If 0 < fc < n, then Ai < 0 < A„. 

Linearized equation and linearized flow 

Linearizing the downward gradient equation (1) at a constant solution u{t) = 
p G Grit/ in direction of a vector field ^ : M —)■ TpM along u leads to the linear 
autonomous ODE 

i{t) = -Am ( 8 ) 

for ^ where A = Ap is the Hessian operator determined by (5). 
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Interpreting —A : TpM —>■ TpM in (8) as a vector field on the manifold 
TpM it is easy to see that the flow generated by —A is the 1-parameter family 
{<i)t}tgR of invertible linear operators on TpM given by 

f, = e-“:=f;fcfli,4r (9) 

If the spectrum of A is given by {Ai,..., A„} C K, then ■ ■ • j C 

(0, oo) is the spectrum of We will use only for 

t > 0. However, the part will be needed for all t € M. Since this part acts 

on the (linear) unstable manifold E~ and is defined by a matrix exponential 
- thus without employing any backward Cauchy problem - there is no conflict 
with Convention 1.5. 

Lemma 2.1 (Linearized flow satisfies linearized equation). Given p G Crit/, 
then the family of linearizations {d(ptip)}t>o coincides with the time-t-maps 
{d>t}t>o associated to the linearized vector field —A = —D{yf)p on TpM. 

Proof. See e.g. [Web06, Le. 2.5]. □ 

Proposition 2.2 (Exponential estimates). Consider the Hessian operator A = 
A~ 0 A~^ assoeiated to a non-degenerate critical point x of Morse index k and 
with induced splitting T^M = E~ 0i?+ provided by (7). Fix a eonstant p in the 
spectral gap [0,(i] of A. In this case there are the estimates 

for every t > 0 and 

II < ||e-‘^” || < < e*'', 

for every t<0. Thus \\e~*^\\ < max{e“‘^'=+C < max{e“*^, when¬ 

ever t > 0 and ||e“*^|| < max{e“‘'^", e”*'*'''} < max{e“*''‘", e*''} whenever t < 0. 

Proof. Eigenvectors associated to Xj and e~*^^ can be chosen equal. □ 

Exponential map 

Given any point q of our Riemannian manifold M, we denote by exp^ : TqM D 
Oi,^ —)■ M the associated exponential map. Here is the open ball whose 
radius tg > 0 is the injectivity radius of exp at the point q. The infimum t > 0 
of Lq over M is called the injectivity radius of M. The map 

exp : TM Z) M 

is called the exponential map of the Riemannian manifold. It is a piece of 
mathematical folklore that one can introduce global maps Ei and E 2 which 
can be viewed as partial derivatives of the exponential map in the direction of 
M and in the fiber direction, respectively. Also derivatives Eij of second and 
higher order can be introduced; see e.g. [Webl4a, §2.1]. 
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Theorem 2.3. Given u € M and ^ G C TuM, then there are linear maps 

: TuM -G Eij(u,^) : T^M x T^M —)■ Tg^p^^M, 

for i,j G {1, 2} such that the following is true. // u : M —)■ M is a smooth curve 
and f,ri are smooth vector fields along u such that ^(s) G for every s, 

then the maps Ei and Eij are characterized^ by the identities 

^exp„(^) = Ei{u,f)dsU +E2iu,f)Vgf 

Vs {Ei{u,f)ri) = Eii{u,f) iri,dsu) + £^i2(m,0 + Ei{u,f)Vsr] 

Vs {E2iu,f)p) = E2iiu,f) iri,dsu) + £^22(m,6 + E2{u,f)VsV- 

These maps satisfy the identities 

Ei{x,0) = E2 {x,0) = 1, Eii{x,0) = E2 i{x,0) = E22{x,0) = 0. (11) 

Tangent coordinates near critical point 

Lemma 2.4. Pick p G Critf and let ip > 0 be its injectivity radius. Then in 
the neighborhood exppOt^ of p in M the solutions to u = —{V f) o u correspond 
via the identity u{t) = expp^(t) precisely to the solutions of the ODE 

i{t)+Aat) = hm), hi0 = -E2ip,0-H^f)e.r>^i + M, (12) 

where f takes values in the open ball C TpM and h{f) abbreviates h{f{t)). 

Proof. Given any smooth curve f in 0^^, set u{t) = expp^(t). Take the deriva¬ 
tive to get that ii = i?2(p, C)? pointwise in t. Thus we obtain that 

E2iP,0i = u = -(V/) o U = -E2ip,0Af - ((V/)exp^« - E2{p,0Af) 
pointwise in t and where the last step is by adding zero. □ 


Lipschitz continuity of the non-linearity 

Lemma 2.5 (Lipschitz near p G Grit/). Recall that f G M) for r > 1. 

Fix p G Grit/ and let ip be its injectivity radius. Set po := tp/4. Then there is a 
continuous and non-decreasing function k : [0, po] [0: oo) vanishing at 0 such 
that the following is true. The non-linearity h : Bpg -G TpM given by (12) is of 
class . R satisfies h{0) = 0 and dh{0) = 0 and the Lipschitz estimate 

\\hiO-hih)\\ <Kip)U-r]\\ 


and its consequence 


^uniquely determined. 


|!d/i(/)u|| < n{p) ||u|| 


(13) 



whenever ||5||, Hrylj < p < po and v G TpM. If f is of class locally near p, 
then there is a constant k* > 0 such that dh satisfies the Lipschitz estimate 

\\dh{^)v - dh{p)v\\ < K* 11 ^ - p\\ ||v|| , 

whenever ||^||, HryH < p < po and v G TpM. 

Proof of the Lipschitz Lemma 2.5. That / G implies h G C"" follows im¬ 

mediately from the definition of h; see (12). Pick G Bpg and set X := p — f. 
It is useful to view f as being fixed and p{X) = ^ + X as depending on X. 
Consider the map defined by 

liX) := h{0 - h{p) = -i? 2 (p,^)-l(V/)exp,C + AC 

+ E2{p,f,+X) ^(V/)expp({-|-X) - A(C + AT) 

and note that 7(0) = 0. Moreover, for any v G TpM it holds that 
dI{X)v = ^\^I{X + Tv) 

= -E2{p,f + X)-^E22{p,f + X) [{VfUp^(^+x),E2{p,f + X)-\ 

+ E2{p ,C + X)-i77(V/)exp^(C+x)i?2(p, c + 
where brackets [...] indicate (multi)linearity. Thus we obtain that 
WHO - h{p)\\ = ||7(X)|| = \\I{X) - 7(0)11 = WdIiaX)X\\ 

< 11-^2 ^llL°°(B3po)ll-®22||L~(B3p)||(V/)expj,||L~(B3p)l|-’^ll 

+ ||£;2-'77(V/)expp7;2 - A|U^(B3p)||7f|| 

=: Hp)W^-vW 

where for instance E 2 abbreviates E 2 {p, •) and L°° denotes the sup-norm. More¬ 
over, existence of some constant a G [0, 1] is asserted by Taylor’s theorem. That 
the function k(p) is non-decreasing is due to the fact that the supremum is 
taken over the closed ball of radius 3p. Indeed the supremum is obviously non¬ 
decreasing if it is taken over larger balls, that is if p grows. For p = 0 all maps 
are evaluated at the origin of TpM, thus «:(0) = 0 since^ 7 ^ 22 ( 7 *, 0) = 0 and 
E2(p,0) = 11 by (11) and since D{Xf)p — A = 0. This proves the Lipschitz 
estimate. Concerning its consequence recall that dh^O^ = limT-_j.o 
and apply the Lipschitz estimate. 

Concerning the Lipschitz estimate for dh note that by (12) we get that 

dh{Ov= ^\^h{f + tv) = E 2 ip, 0 ~^E 22 ip ,0 (V/)expp( 4 ), 7 ^ 2 ( 7 *, 

- 7;2(p,C)-'i7(V/)exp^(5)7^2(p,C)^ + An. 

Since by assumption 71(V/)exp (•) is locally Lipschitz near the origin, so is 
dh. ” □ 

In fact, we don’t even need to use the identity £22 (p, 0) = 0 since anyway (V/)p = 0. 
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Cauchy problem and integral equation 


Proposition 2.6. Given a non-degenerate critical point x € Crit/, let Lx he its 
injectivity radius and consider the ODE (12) in TxM with non-linearity h. Fix 
a radius 

po e (0,miii{l,ia;/4}] (14) 

sujficiently small such that the image exp^ Bp^ C M of the closed ball Bp^ C 
TxM contains no critical point of f other than x. Pick T >0 and assume that 
f : [0,T] —>■ TxM is a map bounded by pg. Then the following are equivalent. 


(a) The map f : [0,T] —> TxM of class is the (unique) solution of 
the Cauchy problem given by the localized downward gradient flow equa¬ 
tion (12) with initial value 'C(O)- 

(b) The map f : [0,T] —>■ TxM is continuous and satisfies the integral equation 

or representation formula 


^{t) = e ‘^7r+C(0)+ / e '^''^TT+h{f{a)) da 
Jo 




(15) 


{t—cr)A 


TT-h{f{a)) da 


for every t G [0, T], In the limit T —>■ oo term three in the sum disappears. 


Proof. Essentially variation of constants; cf. [Tesl2, §9.2] or [Josll, §7.3]. □ 


3 Invariant manifolds 

Suppose X G M is a hyperbolic singularity of the vector field X = —V/, that is 
a non-degenerate critical point of /. The stable manifold of the flow-invariant 
set {x} is defined by 

W^{x) := \q G M \ lim iptq exists and is equal to x}. (16) 

In case of a genuine complete flow one simply considers the limit t —>■ — oo to 
define the unstable manifold W'^{x). Note that - despite the naming - these 
sets are, at this stage, nothing but sets. They are invariant under the flow 
though. A common strategy to endow them with a differentiable structure is 
to represent them locally near x as graphs of differentiable maps and then use 
the flow and flow invariance to transport the resulting coordinate charts to any 
location on the un/stable manifold. To carry out the graph construction one 
introduces in an intermediate step what is called local un/stable manifolds. 

Unstable manifold theorem 

In view of our Convention 1.5 to ignore the backward Cauchy problem, already 
defining W^{x) by the analogue of (16) is not possible. A way out is to consider 
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an asymptotic boundary value problem instead: Consider the set of all forward 
semi-flow lines that emanate at time —oo from the non-degenerate critical point 
X, then evaluate each such solution at time zero and define W'^{x) to be the set 
of all these evaluations. In symbols, the unstable manifold of x is defined by 

lC“(x) := {■u(O) I u : (—oo, 0] —> M, (1), lim u(t)=x}. (17) 

It is non-empty since the constant trajectory u = x contributes the element x. 

Theorem 3.1 (Unstable manifold theorem). Non-degeneracy of x together with 
(f being a gradient flow implies that the unstable manifold W'^{x) is an embedded 
submanifold of M of class C’’ tangent at x to the vector subspace E~ C T^M 
of dimension k = ind(a;) and diffeomorphic to E~. 

Corollary 3.2 (Descending disks). Given x G Crit/ non-degenerate, then there 
is a constant > 0 such that the following is true. 

a) Each descending disk defined by 

Wf{x):=W^{x)n{f>f{x)-e}, ee(0,£,], (18) 

is C’’ diffeomorphic, as a manifold-with-boundary, to the closed unit disk 
jjfe C where k = ind(x). The boundary 

Sf{x) = W-{x)n{f = f{x)-e} 

of a descending disk is called a descending sphere. 

b) Each open neighborhood of x in M, thus each open neighborhood of x in 
lU“(a;), contains a descending disk. 

Suggestion for proof. Apply Theorem 3.1 and the version of the Transversality 
Theorem for manifolds-with-boundary, see [Hir76, Ch. 1 §4], to obtain the C'’ 
manifold structure of the descending disk. Use the Morse Lemma [Hir76], which 
causes loss of regularity, only to control the locus. □ 

The proof of the unstable manifold Theorem 3.1 is a Corollary of the local 
unstable manifold Theorem 3.4 below. The standard argument is to use the 
forward flow to move the coordinate charts provided by Theorem 3.4 near x 
to any point of lU“(a;). This shows that lU“(x) is injectively immersed. Now 
exploit the gradient flow property. To prove Theorem 3.4 we need a backward 
flow, but only on the unstable manifold, which is coherent with Convention 1.5. 

Definition 3.3 (Algebraic^ backward flow on unstable manifold). Given t6(0) € 
W^{x), set q := u{0) and define if-tq '■= u{—t) for t > 0. Note that f/'-tV'-s = 
ijj-t-s and that ip-tq solves the backward time Cauchy problem (1) with initial 
value q. Therefore it doesn’t violate Convention 1.5 if, for any q G lU“(x), we 
use the notation (ptq for any time t, positive or negative. 

While lU“(a;) is obviously backward and forward flow invariant, a descending 
disk still is backward invariant since its boundary lies in a level set. 

® The wording “algebraic” backward flow is only meant to indicate that no backward 
Cauchy problem is involved in its definition. It arises naturally along the unstable manifold 
each of whose points has a past by definition. Thus along VC" (a;) it turns into a genuine flow. 
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Local unstable manifold theorem 


We wish to prove that fF“(a;) carries locally near x the structure of a manifold. 

Theorem 3.4 (Local unstable manifold, Hadamard-Perron [HadOl, Per28]). 
Assuming (H1-H2) in Definition 1.1 there is a constant p = p(A) G (0, such 
that the following is true. A neighborhood ofO in the local unstable manifold 

:= {r?(0) I g : (-oo,0] ^ (12), ^ = o} (19) 


is a graph over the radius R = ^ ball Bj^ C E and this graph is tangent to E' 
at 0. More precisely, there is a C"" map 


= (id, F°°) : B-E-® E+, F^{0) = 0, dF°°(0)=0, (20) 

whose image is a neighborhood of 0 in W^{Q,Bpfi). Moreover, it holds that 
||<(>i''7o|| < for every t <0 and uniformly in 770 € see Figure 5.® 

We sketch the proof of the theorem by the contraction method. Pick an 
element Z- € E~ near the origin. Our object of interest is a backward flow line 
77 : (— 00 ,0] —> Bpg whose value at time zero projects to z_ under 7r_ and which 
emanates from the origin asymptotically at time t = — 00 . For any sufficiently 
small constant p G ( 0 , ^) the complete metric space 


Z- = := {77 G C'°((-oo,0],r,M) 


ll^llexp := supe ||77(t)|| < p| 
^ ^<n J 


carries the contraction given by 


(s) 


z_- 


„-{s-cr)A'' 


TT -h(q(u))da 


+ f e '^'^^TT+h{g{a))da. 
J —00 


( 21 ) 


Its (unique) fixed point is the desired flow line. Define the graph map by 
F°°{z-) := 7r_|_(772_ (0)) as illustrated by Figure 5 and denote by Bf^ C E~ the 
ball of radius R = ^ ^ about 0. 

® Uniform exponential decay: The theorem shows that all backward trajectories which 
remain forever in backward time in a certain neighborhood of the fixed point 0 not only 
converge to 0, but they do so exponentially - even uniformly at the same rate of decay. 
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Stable manifold theorem 

Whereas the stable manifold is easier to define - given only a forward flow - than 
the unstable manifold, the step from local to global is not obvious any more, 
given Convention 1.5. We shall use this oportunity to promote Henry’s [HenSl], 
widely unkown as it seems, argument to pull back the local coordinate charts 
near x in the backward time direction utilizing only the forward flow. 

Theorem 3.5 (Stable manifold theorem). Non-degeneracy of x together with ip 
being a gradient flow implies that the stable manifold W^{x) defined by (16) is 
an embedded submanifold of M” of class C"" tangent at x to the vector subspace 
E~^ C T^M. Thus dim W®(a:) is equal to the Morse co-index n — k of x. 

Remark 3.6. The ascending disk Wf{x) and the ascending sphere S){x) 
are defined as in Corollary 3.2, just replace the superlevel set {/ > f{x) — e} 
by the sublevel set {/ < f{x) +£}. They also have analogous properties, in the 
assertions just replace k hy n — k where n = dimM. 

Theorem 3.7 (Local stable manifold, Hadamard-Perron [HadOl,Per28]). As¬ 
suming (H1-H2) in Definition 1.1 there is a constant p = p(\) G (0, such 
that a neighborhood of 0 in the local stable manifold 

W^{0,Bpg) := Iz G Bpg I ftz G Bpg Vf > 0 and lim ((tz = o| (22) 

L t —^OO J 

is a graph over the radius R = ^ ball B^ C and this graph is tangent to E^ 
at 0. More precisely, there is a C"" map 

= {G°°,id) : B+-^ E-®E+, G“(0) = 0, dG“(0) = 0, (23) 

whose image is a neighborhood of 0 in W®(0, Hp^). Moreover, it holds that 
ll^t^oll < pe~*^ for every forward time t>0 and uniformly in € G°°{B'I^). 

The proof of Theorem 3.7 is by the contraction method. In fact, our 
proof of the Backward A-Lemma, Theorem 1.2, presented below generalizes the 
contraction method proof from invariant manifolds, which is well known, to 
invariant foliations. The details missing in the following sketch of proof can be 
easily recovered by formally setting T = oo in the proof of Theorem 1.2. 

Pick 2 + G E~^ near 0. Our object of interest is a flow line f : [0, oo) —)■ Bpg 
whose initial value projects to z+ under tt+ and which converges to 0 , as t —>■ oo. 
For any sufficiently small constant p G (0, ^) the complete metric space 

Z = Zyp := G GO([0,oo),T,M) | ||e||^,p := supe‘^ \m\\ < p}, (24) 

carries a contraction defined by 

(^'2+C)(i) = + / e~^*~'^^^Tr+h{f{a))da - / 7r-h{f{a))da. 

Jo Jt 
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Moreover, by the representation formula (15) for T = oo the (unique) fixed 
point is our object of interest. For i? = | ^ the map 

: 5+ ^ E-, B+ := {z+ G E+ : ||2+|| < p/2} , 

0+ 7r_ (0)) 

has the properties asserted by Theorem 3.7; see also [CH82, Sec. 3.6]. 

Remark 3.8. Further methods to prove the local (un)stable manifold theorems: 

• Graph transform method: A geometrically appealing method, nicely 
sketched in [PdM82, p.80]; for details see [Shu87]. 

• Irwin’s space of sequences: See former two references. Another excellent 
reference for those who care about details is Zehnder’s recent book [ZehlOj. 

Proof of the stable manifold Theorem 3.5 (Henry [HenSl, Thm. 6.1.9]). 

By Theorem 3.7 the stable manifold is locally near x a submanifold of M 
of dimension n — k where k is the Morse index of x and n = dim M. Thus 
there is a neighborhood E of a; in W^{x) which is represented as a zero set 
{y G U \ h{y) = 0} where U is an open set in M and /i : C/ —>■ is a G’’ map 
such that dhy : TyM —>• is surjective at each point 

y G h-i(o) = c/n E = unw\x). 

Note that here only the first identity is part of the submanifold property of E. 
To obtain the second identity choose U smaller, if necessary, and use the fact 
that a flow line of a gradient flow cannot come back to itself asymptotically. 

Nearby arbitrary elements go of W^{x) one obtains local submanifold charts 
as follows. Pick T > 0 such that (fxqo € ^ and consider the map 

h o ipj' : Ut ——)■ U — )■ Ut := 'Pt 

whose zero set is Ut H W''^{x)\ see Figure 6. The pre-image Ut is an open 
neighborhood of go in M since U is open and pt is continuous. By the regular 
value theorem it remains to show that the map 

d{h o pT)y : TqM Ty=^^gM ^ 

is surjective whenever g G {hoipT)~^{0) = Gt nIF'*(x). Since dhy:,.j.q is surjective 
it suffices to show that d{pT)q is surjective.’^ The following argument avoids 
backward flows. In finite dimensions surjectivity of d{pT)q is equivalent to 
dense range which itself is equivalent, even in the general Banach space case, to 
injectivity of the adjoint (or transposed) operator {dpT\q)*■ But the latter is 
equivalent to backward uniqueness of the {forward Cauchy problem associated 
to the) adjoint equation. In our case A = A* so the adjoint equation is just 
the linearized equation itself. But an ODE associated to a Lipschitz continuous 
vector field exhibits forward and backward uniqueness; see e.g. [AL93]. □ 


^ In infinite dimensions [HenSl] the operator dipTll is not surjective in general, but still 
admits dense image. This suffices to show surjectivity of the composition since k is finite. 
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Figure 6: Henry’s proof: pull back coordinate charts by forward time-T-map 


Convenient coordinates 

Suppose the local setup of (H1-H2) in Definition 1.1. Namely fix a constant A 
in the spectral gap ( 0 , d) and recall that the local flow (ft acting on the ball Bp^ 
in T^M = X = E~ © £'+ is generated by the ODE f + A(f = h{() where the 
non-linearity h is given by (12). Consider the balls B~j^ C E~ and B^ C 
of radius R(x, po,X) = ^ ^ with Bf^ x B^ C Bpo /2 and the C’’ graph maps 

F°° : B]^ —>• E+ and G°° : B^ E~ provided by the local (un)stable manifold 
Theorems 3.4 and 3.7. Use the fact that the graph E°°{Bf^) is tangent to E~ 
at 0 , similarly for Q°°{B^), to see that choosing the radius i? > 0 smaller, if 
necessary, one can arrange that both graphs simultaneously satisfy inclusions 

C {Bf, X B+) , g^{B+) C {Bf, X B+) . 

On the other hand, by Corollary 3.2 and Remark 3.6 there is a (small) constant 
2 ^ > 0 such that there are inclusions of descending and ascending disks 

Wi,cg°^{B+). (25) 

Set 

D- := 7r_TU2“ , D+ := n+W^^, U := D" x D+. (26) 

Then by the graph property of and g°° it holds that E°° {D~) = and 
g°°{D^) = as illustrated by the left part of Figure 7. Consequently is 
a disk, that is a set diffeomorphic to a closed ball. Note that 


w,-, = w-(o,v) = vnw-(o,Bpj, 

W 2 \ = W^(0,V) = vnw^(0,Bpj, 

by negative and positive invariance under (f of and respectively. Fol¬ 
lowing Palis and de Melo [PdM82, Ch. 2 §7] observe that the map 

d(x,y) := (x- G°°(y),z/- E°°(x)) 

defined on Bf^ x B^ satisfies d{0) = 0 and d'd{0) = 1. In particular, it is 
a diffeomorphism locally near the fixed point 0. Choosing R and 2c smaller. 
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D+ 


(V, &*<(>) 





Figure 7: Local (un)stable manifolds get flattened out 


if necessary, we assume without loss of generality that i? is a diffeomorphism 
onto its image. Note that 0) = {x,—F°°{x)) and = {—G°°{y),y). 

Consequently the map d diffeomorphically maps F°°{D~) = to the disk 
D~ and similarly for W 2 ^ and D'^; see Figure 7. Our new local model will 
be the product of disks V := D~ x £>+ acted upon by the new local flow 
(denoted again by (j)t in Remark 3.9) whose local unstable and 
stable manifolds are D~ and respectively. Note that the latter represent 
the descending disk and the ascending disk W 2 ^{x) in this new local 

model. The ODE which generates the flow arises by re-doing Lemma 2.4 
starting from the “Ansatz” u{t) = exp,j, for Obviously the new 

ODE is similar to (12) involving, in addition, the diffeomorphism -d and its 
linearization. Since -d is of class C'*' and defined on the compact set V it admits 
uniform C’'-bounds. Therefore relevant estimates for the new and the old ODE 
are equivalent up to constants. This justifies to simplify notation as follows. 

Remark 3.9 (Coordinates with flattened local manifolds). By the discussion 
above and in order to simplify notation we will, without loss of generality, con¬ 
tinue to work with the ODE (12), that is we omit i? in our formulas, and assume 
that the local (un)stable manifolds and for the (new) local flow (p on 
V are disks D~ C E~ and D'^ C E+. To summarize we assume that 

V := D- X D+ C B]^ X B+C Bi, £>- = D+ = VF^, (28) 

as illustrated by Eigure 1 and where > 0 has been fixed according to (25). 

4 Proof of backward A-Lemma 

Suppose ind(a:) = fc S {1,... ,n — 1}, so the smallest and largest eigenvalue of 
the Hessian operator A = satisfy Ai < 0 < A„. Consider the continuous 
function k{p) with k( 0) = 0 and the Lipschitz constant k* > 0 of the non¬ 
linearity h provided by Lemma 2.5. Pick an exponential decay rate A S (0, d) for 
the elements of the complete metric space to be defined below and consider 
the two constants <5 G (O, min{l, and y G (-^j C (A, d) defined in 

Figure 4. To ensure the second of the two endpoint conditions in (32) set 

Ti =Ti(A,><) :=-i^ >0. (29) 
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Fix p = p{X) e (0,1) sufficiently small such that 

and such that the closed p-neighborhood Up{W^) in X of the descending disk 
is contained in = D~ x D'^ C i?i; see (28). Observe that in this case the 
radius | ball 5+ C E'^ is contained in 0+. Fix T 2 = T 2 (p) > 0 such that 

yggj g^gp 5^ gg^ 

To = To(A,/x(A), k) := max{Ti,T2, 1} > 1. (31) 

Remark 4.1 (Mixed Cauchy problem). The key observation to represent the 
pre-image 4>t~^'E>z_ under the time-T-map (j)T as a graph over the stable tangent 
space T+ is the fact that there is a well posed mixed Cauchy problem. 
Namely, instead of prescribing precisely the endpoint, by the representation 
formula (15) it makes sense to prescribe only the (7r_|_)-part of the initial point, 
but in addition the (7r_)-part of the endpoint. This way one circumvents using 
a general backward flow - only the (by Convention 1.5 legal) algebraic one on 
the unstable manifold is used; cf. [Webl4a, Rmk. 3]. 

We prove below that for each z+ G T+ with ||z+|| < pl2 there is precisely one 
flow line ^ with initial condition 7r_|_^(0) = and endpoint condition 

^(T) e T> 2 _. The latter is equivalent to 

7r_e(T) = z_ A U{T)-z.\\<x. (32) 


The key step to determine the unique semi-flow line ^ associated to the triple 
{T,z-,z+) is to set up a strict contraction 4/^ on a complete metric space Z'^ 
whose (unique) fixed point is More precisely, define 


= zl;- := [i e C%[Q,T],X) : ||?- < p) , 

where z'L '■= (f>-T{z-) and 

llellexp = llellexp.T H - 

and consider the operator defined for ^ G Z'^ by 

(t) ■■= e-*^z+ + e-(‘-‘")^7r+h(^(cr)) da 

+ Z-— f n_h{^{a))da 


(33) 


(34) 


(35) 


for every t G [0,T]. The fixed points of correspond to the desired flow 
trajectories by Proposition 2.6. In Step 1 and Step 2 below we show that 
is a strict contraction on Z^. Hence by the Banach-Cacciopoli Fixed Point 
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Theorem, see e.g. [CH82, §2.2], it admits the unique fixed point and we 

define the map 

-.S^xB+^E-, (0_,z+)^7r_ei:,z+(O)=:GL(^+). (36) 

For latter use we calculate for t € [0, T] and a > 0 the integrals 

T)a 

— ■ (37) 

The proof takes six steps. Fix z_ G S'" and z+ € 5+. Abbreviate 

Step 1. For T > 0 the set equipped with the exp norm metric is a complete 
metric space, any ^ G takes values in V C Bi, and acts on Z"^. 

Proof. For the compact domain [0,T] the space G°([0,T], A) is complete with 
respect to the sup norm, hence with respect to the exp norm as both norms 
are equivalent by compactness of [0,r]. By its definition the subset Z'^ of 
G°([0, T],X) is closed with respect to the exp norm. Since Up{Wf ) C C/p(bF“) C 
V by our choice of p, the elements of Z^ take values in ^ C Bi. 

To see that acts on we need to verify that is continuous and 
satisfies the exponential decay condition whenever f G Z^. By definition 
is a sum of four terms each of which is continuous as a map [0,r] —)■ X. For 
terms one and three continuity, in fact smoothness, follows from definition (9) 
of the exponential as a power series. Continuity of both integral terms, terms 
two and four, uses the same argument. Denote term two by F{t) from now 
on. Continuity of F : [0,T] ^ X and the fact that F{0) = 0 (used in Step 3 
below) both follow from continuity and boundedness of the map h := 7r+ oho^ : 
[0,T] —>• A+ which holds true since ^ : [0,T] —>• A is continuous and bounded 
by definition of and so is the non-linearity h by the Lipschitz Lemma 2.5. 

We prove exponential decay. For t G [0, T] consider the flow trajectory given 
by where := <P-t{z-). By the representation formula (15) it satisfies 


da = 


1 — e 


^ da = -—— 


(ftz'^ 



. (t da 

{t-T)A _ f ^-{t-a)A TY_}i[(j)^z^) da. 


(38) 


Here we used that Z- G Sf = dWf', thus z'^ = (f-rZ-, lies in the (backward 
flow invariant) descending disk IF" C D~. Hence 71+zf = 0 and 77-4)tzT. = 
7r_z_ = 2 _. By (35) and (38) we get for t G [0,T] the estimate 


< ||e“*^z+|| + ||^(?(o’)) - h{(l>az'^)\\ da 

Jo 

+ ^ ||/i(e(^)) - hi<j7,z^)\\ da 


(39) 
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< ||z+|| + n{p) ||^ - da 

+ n{p) lie - ILxp da 

where we used the exponential decay estimates in Proposition 2.2. To the non¬ 
linearity h we applied the Lipschitz Lemma 2.5 to bring in the constant k(p). 
We also used the fact that the norm of a projection is bounded by 1. Moreover, 
we multiplied both integrands by to obtain the exp norms which are 

bounded by p by definition of Inequality three uses that z+ € B'^ and the 
integral estimates in (37). The hnal step is by smallness (30) of p. □ 

Step 2. For T > 0 the map is a contraction on Z^. Each image point 
satisfies the initial condition 7r+ (0) = z+ and for T > Ti also the 

endpoint conditions (32), i.e. it hits 'Dz_ = {z-} x 5+ at time T. 

Proof. Pick ^ 1,^2 € Then similarly to (39) we obtain that 

e‘" II it) - (^^6) it)\\ < f l|/i(ei(a)) - h(6(a))|| da 

Jo 

for every t € [0, T], Use smallness (30) of p and take the maximum over t G [0, T] 
to get that - ^^6llexp < yCi - 6llexp- 

To obtain the identities 7r+ (0) = and 7r_ (T) = z_ just set 

f = 0 in the definition (35) of and use the identities 7r_|_7r_ = 7r_7r_|_ = 0, 
continuity of the exponential series (9), and continuity and boundedness of both 
integrands. Concerning the second endpoint condition in (32) assume T > Ti 
and evaluate estimate (39) at t = T to get that 

II(^'0 (T) - ^-llx ^ 

where the last step is by p < po/2 < 1 and definition (29) of Ti. □ 

Step 3. For T > 0 the map : Sf‘ x 5+ —)■ E~ defined hy (36) is of class G'’ 

and, for each Z- G Sf, the map G^_ := G'^{z-, •) : B'^ -G E~ satisfies 

Gf_(0) = fi-T{z-) =■■ z^, graphGj_ = {5j_..^(0) | 0 + € B+} . 
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Proof. By Step 2 and its proof the map 

xB+xZ^^ Z^, iz.,z+,0 ^ 

is a uniform contraction on with contraction factor (Strictly speaking 
Z'^ depends on 2 _, but the complete metric spaces associated to different Z-^s 
are naturally isomorphic.) Observe that is linear, hence smooth, in 2 ;_ and 
in z+ and of class C"" in f, because h is of class C’’ by the Lipschitz Lemma 2.5. 
Hence by the uniform contraction principle, see e.g. [CH82, §2.2], the map 

e = e'^ : X B+ Z'^ 

f ^ CT (40) 

{z-,z+) ^ 

which assigns to (z_, z+) the unique fixed point of is of class C"" and so 

is its composition with the (linear) evaluation map evo : Z^ —>■ X, ^ i—)■ ^(0), 
and the (linear) projection 7r_ : X —)■ E~. But the composition of these maps 
is by its definition (36). So G^, thus 5, is of class G'’ in z_ and z_|_. 

Consider the flow trajectory r] : [0,T] —>■ X, t >—>■ which runs from 

zf to z_ inside the (backward invariant) descending disk Wf C D~. Hence 
7r_|_?7(0) = 0 and 'K-r]{T) = ry. Thus ry = q by uniqueness of the fixed point. 
Hence G^_ (0) := 7r_^J_ g(0) = 7r_ry(0) = H-zf = zf. To get the desired 
representation of graph Gj_ observe that 


gUz+) ■■= (gJAz+),z+) = (^-4^_,,^(o),^+4^_,,^(o)) = e_,.+(o) (41) 

by Definition (36). The first identity also uses the fixed point property and the 
initial condition from Step 2. The final identity is by 7r_ 0 7r+ = 1.x- D 

Step 4. The map T i—> G{T, Z-, z+) is Lipschitz continuous. If f is of class G^4 
near X, then the derivative T i—)■ -^Q{T, Z-, z+) is also Lipschitz continuous. 


Proof. We prove that Q is Lipschitz continuous in T. Fix T > Tq > 0, z- G Sf, 
and z+ G H+. Consider the fixed point of the strict contraction 

4'^ = defined by (35). The fixed point of 4'’^+'^ is given by 


= e *^z++ / e d '^^"^7r+/i(^’^+'^((T)) dcr 

Jo 




rT-\-T 


(42) 


{a)) da. 


For t € [0,T] and t > 0 we obtain, similarly to (39), the estimate 

|kT+.(^)-^T(^)|| 


lo 

rT 


yj(^T+T(^)) _ ^(^T(^))|| II (gXA (t T)A 


pT pT-\-T 

+ e^^-'^^^\\h{e^^a))-h{eia))\\da + \\hie^^{a))\\ da 
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< k(p) -C^|lco([O.T],X) '")''dcr + ^ e(‘ 

+T 

+ t|Ai| • e^‘“^||z_|| + «;(p) / da 

~<i 

2 1 — 6“"^^ 

< «(p)- - ?^|lco([O.T],X) + ^l^ll + <p) - - - 

- g II? '*' ~? llco([0,T],X) I ■ 

Inequality two uses the Lipschitz Lemma 2.5 for / and the exponential estimates 
in Proposition 2.2. Moreover, we used the fact that < 1 to obtain 


1 - 


= e ds < T. 


The identity 


^tA 


-1 = A- 


^aA 


da 


(43) 


(44) 


follows by definition (9) of the exponential as a series. Together with the expo¬ 
nential estimates in Proposition 2.2 we get that 


-11 


A- 


^(tA 


da 


< |Ai| / e da < |Ai| r. 


(45) 


Coming back to inequality two, it remains to explain the estimate for term four. 
Here we used that g takes values in Bi by Step 1. 

Inequality three uses (37) for the integrals and the fact that || 2 _ || < 1. Inequality 
four uses estimate (43) and the smallness assumption (30) on p by which k(p) < 
I. Now take the supremum over t g [0, T] to obtain that 

II? ~? Ilco([o.r],jf) — 


with constant ci = 2(|Ai| -I- 1). By (41) this shows that 


- gl_ (z+) = (0) - ^^(0) < ciT 


(47) 


for all T >To and r > 0. In other words, the map T ^ Q{T, z-,z+) = g^_ {z+) 
is Lipschitz. The difference is illustrated by [Webl4a, Fig. 5]. 

IFe prove that the map T i—>■ -^Q{T, Z-, z+) is Lipschitz continuous. Set 

0^(t) = 0f_ (t) := £ 4^+1^ (t) (48) 

for every t g [0,T]. Calculation shows that this derivative is given by 


0"^(t) = /" e o 0^(cr)) da-I-H e 2 :_ 

Jo 

- e"(‘"^)'^"7r_/i(^^(r)) - ^ (d/i|^T(^) o 0 ^(ct)) da. 
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Since ^Q{T, Z-, z+) = 0^(0) by (41), it remains to show that the map T i-> 
0^(0) € X is Lipschitz continuous. By definition of 0^ we get the identity 

— 0^(t) = f o (0 ^+'^(ct) — Q^{a)) da 

Jo 


+ / e (c?/i|jT+x(^) — c?/i|jt(ct)) o 0’^(cr) da 

Jo 

+ {e^^~ - z- 

- {e^^~ - l)e-(‘-^)^"7r_/i(^'^+^(T + T)) 

- {ne^^T + r)) - fe(g^(r))) 

- 7r_d/i|jT+T(^) o (0’^''"'^(cr) — 0 "^((t)) da 




(t a)A 7r_ — dh|jT(cr)) o 0 ^((t) dcr 

g-(*-<^)^ 7r_d/i|^T+T(^) o Q^~^^{a) da 


for all t G [0, T] and t > 0. To obtain lines one and two add zero; similarly for 
lines four and five and lines six and seven. Abbreviate the T],X) norm 

by IMIc“ • Combine line one with line six and line two with line seven to obtain, 
similarly as above and again abbreviating ci = 2(|Ai| + 1), the estimate 


< k{p) || 0 ^+" - 0 ^ 11^0 ^ ^ 

+ - ellco l|0nico da + £ da^ 

+ t|Ai|^ ||z_|| +r |Ai| k{p) + t)|| + njp) ||g^+^(r + r) -g^(r)|| 

pT +T 

< k{p) ||0"^~''’' — 0 ^||(^o — I" K *( cir)ci — h T |Ai| (|Ai| + 1) 

'"-'T p p 

1 - 

+ k.{p)t (ci + 2A„ + 1 ) + k{p)ci - 

d 

- 3 ||0^’^^ - 0^|lc« +'^ ( ^ ^ + -^ + (ci+2A„ + l) + ci^ 

for all t G [0,T] and r > 0. Inequality one uses the exponential estimates in 
Proposition 2.2 and the Lipschitz Lemma 2.5 for dh. To obtain line three we 
used (45) and we added zero in the form of li(0), thus bringing in k{p). 
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Inequality two uses the following arguments. Estimate the first pair of integrals 
by I using (37); similarly for the second pair. Apply estimate (46) and recall 
that 0 - G Bi by our local setup. Use again (46) to conclude that 

lle^WlI = IIb'G«)|| = lim „ ,49) 

whenever t G [0,T]. Thus ||0^||cg, ^ ci and, of course, the same is true when 
T is replaced by T + r. By Step 1 the elements of the complete metric spaces 
Z'^ (and 2'^+'^) take values in Bi. We also used that k{p) < 1 by (30) and that 
< I 2 ^/jg estimate for the difference ^^+'^(T + r) — ^'^(T) G X in line 
four will be carried out separately below; see (51) for the result. 

Inequality three uses the smallness assumptions (30) on p and estimate (43). 
Now take the supremum over t G [0,T] to get the estimate 

II® llco([0,T],JC) — 

for all T >Tq and t > 0 and where C = (^(Ai, A„, k*, /i“^) is a constant. So 

II ^g{T + T, z_, z+) - ^g{T, z+)II = ||0^+"(O) - 0^(0) II < Cr 

which means that the map T i-G -^g(T, z-, z+) is Lipschitz continuous. 

As mentioned above it remains to estimate the norm of the difference: 

+ r) - f{T) = (e-"^ - ll) e-^^z+ 

+ [ (/i(^^+^(ct)) -/i(C^(cr))) dcr 

+ f — ]l) dfj 

Jo 

pT +T 

+ J e-^^+^-'^^^Tr+h{f+^a)) da. 

Note that we added zero to obtain terms II and III in this sum I+II+III+IV of 
four. We proceed by estimating each of the four terms individually. 

I) Concerning term one the identity for A~^ corresponding to (44) shows that 

||(e-^-ll)e-^^z+|| = r -A+e-^^+^^^^z+ds 

Jo 

< ||A+|| • II2+II • [ e"(^+'^^'"(is 

Jo 

< Xne~'^^- —-< T\ne~'^^ 

h 

where the last step uses (43) and A„ is the largest eigenvalue of A; see (6). It 
might be interesting to see how the in infinite dimensions infinite quantity ||A“'"|| 
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can be avoided, also in III) below; cf. [Webl4a, p. 954]. 
II) For term two we get the estimate 



,-(T+r-<T)A 




da 


< ciTK(p)e-(^+")'^ [ e’^^^da 

Jo 

where we used (46) and smallness (30) of p. 

Ill) Term three requires similar techniques as term one and we obtain that 


r (e--^ - 11) e-(^-‘^)^^+/i(C^(a)) 

Jo 

< [ [ A+e"(^+^"'^^^’^7r+/i(^^(cr)) 

Jo Jo 

f e-('^+"-‘")^dsdCT 


da 

ds da 


< A„k(p) 


/o ^0 


1 _ p-TfJ. 1 _ p-Tfj, \ 

\ / \ ^ ^ 1C An 

= Xni^ip) -< —r 


fl /i 5 

where we used once more that ||^^(ct)||x < 1 by Step 1. We also applied (43). 
IV) For term four we get the estimate 



e-('^+^-'")^7r+/i(e^+^((T)) 


1 - 

da < k{p) - 


< r 


where we used once more the estimates K{p) < 1 and (43). 
To summarize, the above estimates show that 


e^^T + T)-e{T) <r(ci+2A„ + l) 


for every r > 0. This concludes the proof of (50) and therefore of Step 4. 


(51) 

□ 


Step 5. (Exponential convergence) Given 2 _ € S'" and € i?+, then it 
holds that \\Q°°{z+) — 5j_(z+)|| < e~'^s for every T > T 2 . 


Proof. Assumption T > T 2 will be used in (52). Fix z+ € and z_ € S“ 
and consider the fixed point of 'I'^_ on Z'^ and the fixed point 

^ of 4'^^ on Z-, cf. proof of Theorem 3.7. Since Q'^_{z+) = 'C^(O) by (41) 

and Q°°{z.^.) = ^(0) it remains to estimate the difference ^(0) —^^(0). Motivated 
by [Webl4a, Fig. 6] the key idea is to suitably decompose the interval [0,T]. 
The decomposition is illustrated by [Web 14a, Fig. 5] which is surrounded by 
ample explication. Set 

||C-^^|lco := sup ||C(i)-?^(i)|| ■ 

te[0,T/2] 
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Pick t € [0, use the representation formulae for ^ and the Lipschitz 
Lemma 2.5, and add several times zero to get that 


< 


/ g (t o’Ia' ||/j o — h o ^^(cr)|| da 

Jo 

— t\ 

+ J +J I ||/l o ^(cr) — h o ^’^(cr)|| d(T 





Z- 


^{t — cr)A 


||ho,f(CT)|| da 


< + k{p) 








e(t (II^^(ct) - da 


K(rt ||«’'W|| <i<r + „(rt||{||„^yj 




The domain of integration /^2 ^ misprint. 

To continue the estimate consider the last three lines. Now we explain 
how to get to the corresponding three lines in (53) below. Concerning line 
one note that since t € [0,T/2]. Now use the choice of T 2 

and for the two integrals use (37). In line two use that < 1, that 

ll^^(cr) — (t)a{z'^)\\ < pe~'^^ by definition of and that 


j-t{z-) da = 


3T 
'* 8 , 


l>_f_T{z-) dt < 


— t\ 


dt. 


(52) 


Here the identity is by the change of variables t := —a+^T. Set := To 
see the inequality consider the trajectory p{t) := 4>r-T/s.{z-) = 4 >t{z-) defined 
for r € (—00,0]. Observe that p{t) —>• 0, as r —>■ — 00 , because 2 _ lies in 
the descending sphere S'” = SIT” by assumption, thus lies in the backward 
flow invariant set ^-r/glT”. Consequently the whole image of p is contained 
in (j)_T/%W^ C IT” C D~ C E~. Note that 7r_r7(0) = zl and consider the 
unique fixed point of the contraction $ 2 * on Z~ defined by (21). Backward 
uniqueness then implies that p = pz*_. Thus p € Z~, so ||?7(t)|1 < pe^^ < e^^ 
for every t < 0 and this proves the inequality (52). To summarize, line two is 
bounded from above by 


k{p) ■ p 


^ — <7\ 


da - 


3T 
' 8 




dt 


< k{p) • ^ (e ^2 -I- e ^8^ 


Concerning line three note that ||C^(cr)|| < 1 since the elements of Z'^ take 
values in Bi and that ||CI|exp < P < 1 by definition (24) of Z. Now calculate 
both integrals (in the second one use < 1). 
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Overall we get that 


2k(p) rpX 

+ — 

A 

T''T X 

+ -^e~^ 4 + -^e~^ 2 

/i A 


(53) 


for every t G [0,T/2]. Here the last inequality uses the estimate ^ = A + <5 > A 
and smallness (30) of p. Take the supremum over t G [0, T/2] to obtain that 

||C-Calico < (54) 

II "^T/2 7 

Since \\g°°{z+) - gf_{z+)\\ = ||^(0) - C^(0)|| < this proves Step 5. 

□ 


Step 6. (Exponential convergence) Pick 2 _ G 5“, z+ G i?+, v G E~^, then 


dgj_{z+)v - dg°°{z+)i 


< c*e 




for some constant Ci,{k,i,, 1/5, 1/A) and every T >Tq whenever / is of class 
near x. 


Proof. One would expect that the proof of convergence of the linearized graph 
maps should use convergence of the graph maps themselves. Indeed esti¬ 
mate (54) above is a key ingredient. Fix T > Tq, Z- G Sf, 2 + G i?+, and 
V G E+. We proceed in three steps. Step I and Step II are preliminary. 

I. For r > 0 small consider the unique fixed point := € Z'^ of 

the contraction defined by (35). Set := f,Q = Cj_,z+- By (35) 


CrW = e *^(z+-bTv)-b / e ‘^)^7r+/i(C^(cr))dcr 

Jo 

+ z-— 7r_/i(^J((j)) dcr 


(55) 


for every t G [0, T], By the proof of Step 3 the composition of maps 

^l-,Z++TV Cj_,Z+-|-TJ)(0 


is of class C^. Therefore the linearization 


x^{t) = xl 


2+; 


/ N d 

Tr 


(L 


T — 0 


,2++ri? 


it) 


(56) 
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is well defined. For every t G [0,T] it satisfies the integral equation 


(t) = e + [ e '^'>^TT+ (dh| 5 T(„) o Xj( ct)) da 
Jo 

- 7r_ (d/i|^T(^) oxj(cr)) dcr 

and therefore by straighforward calculation the estimate 

e‘^||xj(t)|| < < 2||r;|| 

for every t G [0,T]. Use (41) to see that 

xj(0) := £\^^o^L,z++rviO) =dGLiz+)v. 


(57) 


(58) 


(59) 


II. Concerning the local stable manifold observe the following. For r > 0 
small consider the unique hxed point ^z^+tv of 4'z++ru on Z defined by (24). 
It satisfies the integral equation 


^z++rv{t) = e *^(z+ + ru)+ fed '^'>'^7r+h{^z++rvicr)) da 

Jo 




Tr_h{^z++Tvicr)) da 


(60) 


for every t > 0. Hence the linearization 


Xy(t) = Xz:^.y{t) := — 


^Z^+Tv{t) 


T — 0 


satisfies for every t > 0 the integral equation 

Xy{t)=e~^^v + J e~d-'^^^Tr+o Xy{a)^ da 

- (^T:_dh\^^^(a )0 Xy{a)^ da 

and therefore by straightforward calculation the estimate 


e^^\\Xy{t)\\ < ||X„ 


lexp < 2||^^|| 


for every t > 0. Since Q°°{z+) = ^z+(0) we get that 


(61) 


(62) 


^-(0) := £ L=o ^-++-(0) = dg^{z+)v. (63) 

III. Abbreviate := and ^ := To estimate the difference 

Xy — Xy consider the corresponding integral equations (57) and (61). Add zero 
and apply the Lipschitz Lemma 2.5 for dh to obtain that 
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||d/i|^T(^)Xj(cr) - dh\^(^^)Xy{a)\\ 

= II {dh\iT^^) - dh\^(^)) Xj (a) + dh\^t^^) (xj (cr) - Xy{a)) || (64) 

< K* ||C^(ct) -^(ct)|| • ||Xj(cr)|| +k{p) ||Xj(cr) -X„(cr)|| 


for every a € [0,T]. Here we applied estimate (13) to get that 

||(ih| 5 (^) o X„(ct)|| < k{p) ||X„(cr)|| . (65) 

Pick t € [0, ’^] and apply estimates (54), (58), and (62) to obtain 

||^J(^)-^^ 

< 


e ^)^(^/t«J|g^(a) -^(cr)|| ||xj(cr)|| +k{p) \\X^ (a) - Xyia)\\^da 


<|e-TX/s <2e-''^|| 


+ e(* '")''(K»||^^(CT)-C(cr)|| ||xj(cr)|| +K(p)||Xj(cr)-X^(cr)||)(icr 

pT poo 

+ ||xj(cr)|| d(T + K(p) Il-’^«(f7-)|i da 

^ ^ <2e-‘^^||i;|| 

f da + / da) 

0 Jt I 


bt _^ A 11 

< -K*e 8 ||w 


+ K(p)||xJ-r„||^„^ (y e-da+ y e(‘-")''da 


pOC) 

+ 4k(p)||v|| / da 


< K 


g iA g (A 


d A 

_A 1 


4k(p) ^ 


+ k(p) ( ^ + )^ ) II^J “ 






\xT - X„ 


where c* := 2 k*(| + A) + 1/4. In inequality three we calculated the integrals 
and estimated 4 < A. The final inequality is by smallness (30) of p. Now take 
the supremum over t G [0,T/2] and use the resulting estimate to continue 


dGL{z+)v 


dg°^{z+)v\\ = ||xj(0)-x, 

< 


- 


xj - X, 


r'O 

'-'T/2 


This concludes the proof of Step 6 and Theorem 1.2. 


□ 
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5 Invariant stable foliations and induced flow 


Given x G Grit/ non-degenerate, set c := f{x). Gonsider the local model on 
V = W 2 ^ X W 2 t; C TxM provided by (H3) in Definition 1.1 and recall that the 
assertions of the backward A-Lemma, Theorem 1.2, hold true for each e G (0, <;). 

Invariant stable foliations 

Theorem 5.1 (Invariant stable foliation). Given x G Grit/ non-degenerate, set 
k = ind(a;) and c = f(x). Then for every sufficiently small e > 0 the following 
is true. Pick a tubular neighborhood 'D{x) (associated to a radius x normal disk 
bundle) over the descending sphere Sf{x) in the level hypersurface {/ = c — e}. 
Denote the fiber over q G Sf{x) by Dq{x); see Figure 3. Then the following 
holds for every suffciently large t > 1. 

(Foliation) The subset Nf of M is compact and contains no critical points 
except X. Moreover, it carries the structure of a continuous codimension 
k foliation^ whose leaves are parametrized by the k disk ip-rWf{x) = 
Nf n lT“(a;). The leaf Nf,{x) over x is the ascending disk Wf{x). The 
other leaves are the codimension k disks given by 

= {(pT~^Dq{x) n{f <c + , q'^ := (p^Td, 

whenever T > t and q G Sf{x). 

(Invariance) Leaves and semi-flow are compatible in the sense that 
p&Nfiq'^) =P G VcrG[0,r-T). 

(Contraction onto ascending disk) The leaves converge uniformly to the 
ascending disk in the sense that 

dist {Nfiq'^), Wf{x)) < (66) 

for all T > T and q G Sf{x).^ If U is a neighborhood of Wf{x) in M, 
then A(|(e,T*) C U for some constant t*. 

(Shrink to critical point) Assume U is a neighborhood of x in M. Then 
there are constants £* and r* such that C U. 

Theorem 5.1 is illustrated by Figure 8. Theorem 1.2 asserts that some neigh¬ 
borhood of X is the union of images of (n—fc)-disks. The first step (Gorollary 5.3) 
is to show that these images are disjoint and compatible with the flow. The ac¬ 
tual proof of Theorem 5.1 coincides word by word with the infinite dimensional 
version [Webl4b, Thm. G]. To see this coincidence set i? = | and impose on 
£ > 0 the condition Wf C Note that C D+. 

® For the precise degrees of smoothness away from Wf{x) see Theorem 1.2. 

® Here dist denotes the infimum over all Riemannian distances of any two points and the 
constant A £ (0,d) has been fixed in (H2). 
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Local non-intrinsic foliation 

Definition 5.2. A set T) has the no return property with respect to a forward 
flow (/) if I? n = 0: equivalently n I? = 0, whenever t > 0. 

Clearly any subset of a level set has the no return property with respect to 
a gradient semi-flow. Recall that denotes the stable manifold graph map 
provided by Theorem 3.7 and Q is the graph map which appears in the backward 
A-Lemma, Theorem 1.2. 

Corollary 5.3 (to the Backward A-Lemma). Under the assumptions of Theo¬ 
rem 1.2 and with the additional assumption that has the no return prop¬ 

erty the following is true. The subset 

F = := (img;uim^;“) C D- x B+ C V 

of X = TxM carries the structure of a continuous foliation of codimension k = 
ind(a;).^° The leaf over the singularity 0 is given by the subset F(0) := 
of the ascending disk = D'^ and the leaf over the points 

:= f-ra = QI (0), T > Tq, a € Sf, 

are given by the graphs F{a^) := Qff{B^). Leaves and flow are compatible in 
the sense that 

z e F{a^) ^ (j)^z G F{4)^a^) 

whenever the flow trajectory from z to 4>„z remains inside F. 

Proof of Corollary 5.3. By Theorem 1.2 it remains to check that the sets F{a^) 
and F{fi^') are disjoint whenever (3^, that is {T,a) (S',/3)- Assume by 

contradiction that Qf^{z+) = Q^{z+) =: z for some z+ € B~^. Then by (41) the 
point z is the initial value of a flow trajectory ending at time T on the fiber 
Da and also of a flow trajectory ending at time S on Dp. By uniqueness 
of the solution to the Cauchy problem with initial value z the two trajectories 
coincide until time min{T, S}. If T = S, then a = ft and we are done. Now 
assume without loss of generality that T < S, otherwise rename. Hence 
meets T>ct at time T and T>p at the later time S. But this contradicts the no 
return property of D. 

We prove compatibility of leaves and flow. The fixed point 0 is flow invariant. 
Its neighborhood A(0) in the ascending disk is trivially flow invariant in the 
required sense, namely up to leaving F(0). Pick z G F{a^) := Q'fl{B^). By (41) 
the point z is the initial value of a flow trajectory ending at time T on the fiber 
Va. Assume the image ^^([0,r]) = is contained in T’ := im^ U 

and pick cr G [0,T— Tq]. This implies that z+ := G H+. The flow 

line (l)[o,T-iT]4'aZ runs from f^z to ipTZ G Da- Thus by uniqueness this flow 

For the definition of foliation see e.g. [Law74] or [Pes04, Sec. 4.1] and for the precise 
degrees of smoothness away from W^{x) see Theorem 1.2. 


30 



line coincides with the fixed point ^a,z1 of fh® strict contraction But 

= Ca.2^(0) is equal to Q'^~’^{z+) again by (41). Thus 

<kaz e Ql-%B+) =: F{a^-) = F{cj)„a^) 
by definition of F and a^~'^. □ 


Induced flow — Dynamical thickening 

Theorem 5.4 (Strong deformation retract). Consider a pair of spaces {Nf, L|) 
as defined by (3-4)- Then the following is true for all pair parameters £ > 0 suf¬ 
ficiently small and r > 1 sufficiently large. Firstly, the pair strongly deformation 
retracts onto its part in the unstable manifold, that is onto 

{D,A) = {Dl,A%) := {Nf,CW-{x),L%CW-{x)). 

Moreover, this pair consists of the closed disk D = (p-rWfi{x) whose dimension 
k is the Morse index of x and an annulus A which arises by removing from D 
the smaller open disk vatip_ 2 TWf{x); see Figure 3. 

Corollary 5.5. For any pair of spaces in Theorem S.f it holds that 




Z ,£ = k := ind(x), 
0 , otherwise. 


(67) 


To construct a deformation to prove Theorem 5.4 there is the immediate 
temptation to use the already present forward flow (pf. Unfortunately and ob¬ 
viously, this only works along the stable manifold where indeed the flow moves 
any point into the unstable manifold, as time t —)■ oo. Along the complement 
of the stable manifold this does not work at all.^^ However, given the foliation 
in Theorem 5.1 of Nf in terms of graphs, it is then a natural idea to turn this 
foliation into a dynamical foliation in the sense that each leaf will be endowed 
with its own flow. The natural candidate for these flows is our given flow pt on 
the ascending disk Wf{x) = N^r\W^{x) transported to any leaf via conjugation 
by the corresponding graph map; see Figures 8^^ and 9. This way the set Nf 
turns into a disjoint union of copies of the dynamical system {Wf{x),p). So 
one could call this procedure a dynamical thickening of the stable manifold; 
see [Web] for an application to a classical theorem. 

Proof of Theorem 5.4- Throughout we work in the local model provided by Def¬ 
inition 1.1. In particular, we will use these notation conventions in our con¬ 
struction of a strong deformation retraction 6* of IV onto its part D in the un¬ 
stable manifold. As pointed out above on the stable manifold the forward flow 

This is known as discontinuity of the flow end point map on unstable manifolds and 
obstructs simple geometric proofs of a number of desirable results such as the one that the 
unstable manifolds of a Morse-Smale gradient flow on a closed manifold naturally provide a 
CW decomposition; see e.g. [Bot88,BH04,Nidi,BFKll,Qinll]. 

In Figure 8 we use simultaneously global and local coordinate notation for illustration. 
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Figure 8: Dynamical thickening {N,9) = of ascending disk W^{x) 


{0i}tG[o,oo] itself does the job. Indeed (j)oo pushes the whole leaf iV(0), that is the 
ascending disk by Theorem 5.1, into the origin - which lies in the unstable 
manifold. Since (jit restricted to the origin is the identity, the origin is a strong 
deformation retract of N{0). If the Morse index k is zero, then N = N{0) and 
we are done; similarly for k = n. Assume from now on 1 < fc < n — 1. The main 
idea is to use the graph maps Q°° and Gj provided by Theorems 3.7 and 1.2, 
respectively, and their left inverse 7r+ to extend the good retraction properties 
of (pt on the ascending disk N{0) to all the other leaves N{a'^) provided by 
Theorem 5.1. 

Definition 5.6 (Induced semi-flow - Dynamical thickening). By Theorem 5.1 
each z G N = N^’'^ lies on a leaf N^a^) for some T > t and some a in 
the ascending disk S'" and where := (p-xOi. Then the continuous map 
0 : [0, oo] X —>• A defined by 

9tZ := o TT+ o pt o o 7r+(z) = g'^ o o 7r+(z), (68) 

is called the induced semi-flow on TV; see Figure 9. It is of class away 
from the stable manifold. Here (68) simplifies as g°° = {0,id) and tt+\e+ = I- 

Observe that while 9 apriori only takes values in the image F D TV of the 
graph maps, it does preserve the leaves of F; see Corollary 5.3. Now continuity 
on [0,oo) X TV follows from continuity of the maps involved. Given z G TV(a^), 
set := TT+z. Since z+ lies in the stable manifold by construction of our ’’flat” 
coordinates, we obtain that z+(t) := 4>tz+ —)■ 0, as t —)■ oo. So firstly the limit 

9^z := lim 9tZ = g^ lim z+(t) = t/J(0) = 

t—¥00 t—yoo 

exists and lies in the unstable manifold indeed. Here we used continuity of g^ 
and the flnal identity holds by Theorem 1.2. Secondly 9tz —>■ 9aoZ, as T —>■ oo. 
This shows, firstly, that the map 

9^-.N^D, D = P-rW^ ^ {0} U ([t, oo) x 5“), (69) 
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Figure 9: Leaf preserving semi-flow 9t 


is a retraction and, secondly, that the map 9 extends continuously to [0, oo] x N. 
The fact that the origin is a fixed point of 7r-(_ and (j)t implies that 

9ta^ = Ga o(/<t O7r+(0) = 0 q(O) = . 

Hence 9t\D = ido, for every t € [0,oo]. It remains to show that 9t actually 
preserves N. It suffices to show that 9t preserves each leaf of the foliation 

iV = iV(0)U IJ N{a^). 

T>t 

OL^S'^ 

In contrast to the infinite dimensional case [Webl4b] a simple compactness 
argument will do. Note that any leaf, other than A^(0) = is of the form 

N{a^) = n {/ < c + 4, dN{a^) = 0^(5+) n {/ = c + e}, 

whereas iV(0) = G°°{B^) n {/ < c -I- e} = and dN{Q) = S^. Now since the 
boundary of iV(a^) lies in a level set of / and 9s preserves the graph of we 
only need to show that there is a constant —/r < 0 such that 

f{dtz) < -M < 0, Vz e dN{a^), Va^ G D. (70) 

t=o 

This means that the 9 flow is inward pointing along the boundary of each leaf 
and then we are done. But (70) holds true for some constant, call it —2/i, along 
the (compact) boundary of the leaf N{0), just because 9t = 4 on A^(0) and 
/ strictly decreases along its downward gradient flow - unless there is a critical 
point which it isn’t on S§. Compactness of the leaf space D and of the boundary 
of each leaf, together with continuity of the maps whose composition is 
then implies that (70) holds true for all nearby leaves. To restrict to nearby 
leaves just fix e > 0 sufficiently small and t > 1 sufficiently large. □ 


d 

dt 
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Global foliation 


To extend the foliation of the neighborhood of x all along the stable manifold 
W^{x), that is in the backward time direction, according to Convention 1.5, that 
is without using the backward flow, is an interesting open problem - even more 
so in heat flow situations such as [Webl4b]. 
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